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Abstract: 

We introduce a new method to determine the redshift of unknown-redshift BL Lac Objects. The method relies on si- 
multaneous multi-wavelength (MWL) observations of BL Lac objects in optical, X-ray, HE (E^lOO MeV) gamma-rays 
and VHE (E^lOO GeV) gamma-rays. It involves best-fitting spectral energy distribution (SED) from optical through HE 
gamma-rays with a Synchrotron-Self-Compton (SSC) model. We extrapolate such best fitting model into VHE regime, 
and assume that it represents the intrinsic emission of the object. We then compare the observed VHE flux which has been 
affected by the interaction with Extragalactic Background Light (EBL). Constraining the measured vs intrinsic emission 
leads to the determination of gamma-gamma opacity. Comparing the obtained opacity with the predicted opacity based 
on EBL model, we obtain the redshift of the photon source. 

Keywords: BL Lacertae objects: general - diffuse radiation - gamma rays: galaxies - infrared: diffuse background 
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1 Introduction 

In a recent paper 0, some of us (NM, MP and VT) have 
shown how it is possible to use simultaneous multiwave- 
length (i.e. optical, X-ray, HE 7, and VHE 7) data to ob- 
tain a measure the Extragalactic Background Light (EBL). 
EBL is the integrated light produced and re-processed in 
all the history of the universe. Knowledge of the EBL is 
particularly important in studying the VHE 7-ray emission 
from extragalactic sources (typically, blazars), because in- 
teraction of TeV photons with the EBL attenuates their ob- 
served intensity while the resulting cascade will enhance 
the observed GeV emission. These effects are clearly more 
marked for more distant sources, because of the larger op- 
tical depth, t 77 , for the 7vhe — 7ebl interaction: for a 
given EBL model, the quantitative details of the modifica- 
tion to the VHE 7-ray spectrum emitted by a source will 
only depend on the redshift of the source. 

By reversing the line of reasoning in Q, if we assume to 
have a reliable EBL model, then simultaneous multi-A ob- 
servations of a source may allow for an independent deriva- 
tion of the source redshift, z. For instance: (1) low- and 
high-energy (eV, keV, GeV) spectral-energy-distribution 
(SED) data (that are unaffected by EBL) are fitted with an 
emission model; (2) the model SED is extrapolated into the 
TeV domain, and the extrapolation is assumed to represent 



the insinsic (emitted) TeV emission; (3) the 'emitted' VHE 
spectrum is compared with the measured VHE spectrum, 
and the discrepancy between the two is taken to be a mea- 
sure of t 77 ; (4) given a specific EBL model, r 77 can be 
solved as a function of the source distance (i.e., for a given 
cosmology, of z). 

In what follows we present preliminary results of an ap- 
proach to the problem of distance determination for a 
blazar, which is slightly different from the one above and 
in which we employ a dynamical EBL correction during a 
SED fit in which z also appears as a parameter. Just below, 
we describe the source that we are using for our analysis. 
Then we provide some details about the numerical algo- 
rithm that, elaborating what is done in J4), is used to per- 
form the numerical fit to a chosen SED model which will 
return, inter alia, an estimation of the source redshift. 

2 The source and the data 

In this preliminary study we choose a source whose red- 
shift is actually known. By independently determining z 
with the method detailed below, and comparing the result 
with the known value, we will gather an indication on the 
reliability of our method. The source we consider here is 
PKS 2155-304, whose redshift is z = 0.12. We apply our 
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Figure 1: The dataset of PKS 2155-304 fl]. 



method to the simultaneous data presented in fl] and shown 
in fig. 1 . 



3 The method 



Figure 2: The flow chart of the code 

ergy if the z that is currently used in the \ 2 -minimization 
procedure falls within to tabulated values. The flow-chart 
of the code that we employed is shown in fig. 2. 



We assume that the description of the emission of 
PKS 2155-304 can be understood within a one zone 
Synchrotron-Self-Compton (SSC) model: the two peaks, 
synchrotron and Compton, characteristic of these mod- 
els are here understood within the framework detailed in 
0. Although different models have also been employed 
to describe the emission of PKS 2155-304, as discussed 
in (21 the model that we are choosing can also be used 
to fit the data and does not provide quantitatively signifi- 
cant discrepancies in the characteristic parameters that de- 
scribe the source emission; in particular the fit of the data 
to the [3] faithfull SED can be done as described in ||4] 
using the Levenberg-Marquardt method 0, which is an 
efficient standard for non-linear least-squares minimiza- 
tion that smoothly interpolates between the inverse Hessian 
method and the steepest descent method. 

In [4] the known redshift of the source was used to correct 
VHE data, before performing the fit. In the present analy- 
sis, instead, z is what we wish to determine. We can obtain 
our result by a refinement of the fit algorithm described in 
||4), in which z is added as a parameter to be estimated by 
the fit procedure. Of course, while the \ 2 minimization 
procedures explores the parameter space, different z values 
will also be explored. At each step in which z changes, we 
need to apply an EBL correction appropriate for the z that 
is being considered. In other words, EBL correction has to 
be performed dynamically during the fit. The EBL model 
we have chosen is the one by |6|, which represents the min- 
imum EBL level because it is based on measured galaxy 
counts. EBL corrections have been added to the code for 
several possible values of z. Linear interpolation is used to 
determine the value of the EBL correction at a given en- 



4 Preliminary results 

Our results show that the procedure described above is able 
to fit the SED parameters and the redshift of PKS 2155- 
304, and that all the derived values are compatible with the 
published values for the source parameters. 

The statistical analysis of the significance of the fit is cur- 
rently underway. Preliminary results suggest that blazar 
SED fitting to simultaneous multi-A data with dynamical 
EBL correction can indeed be used to provide redshift esti- 
mations for these sources. A more in-depth analysis will be 
presented elsewhere J7J. It will also be important to check 
the stability of the results versus the choice of both the 
emission model and the EBL model - especially in cases 
where no other methods for redshift estimation can be used. 
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